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ABSTRACT 

We present a catalog of 8358 sources extracted from images produced by the Bolocam Galactic Plane Sur- 
vey (BGPS). The BGPS is a survey of the millimeter dust continuum emission from the northern Galactic 
plane. The catalog sources are extracted using a custom algorithm, Bolocat, which was designed specifically to 
identify and characterize objects in the large-area maps generated from the Bolocam instrument. The catalog 
products are designed to facilitate follow-up observations of these relatively unstudied objects. The catalog 
is 98% complete from 0.4 Jy to 60 Jy over all object sizes for which the survey is sensitive (< 3.5')- We 
find that the sources extracted can best be described as molecular clumps - large dense regions in molecu- 
lar clouds linked to cluster formation. We find the flux density distribution of sources follows a power law 
with dN/dS oc ^-^ '^i" ! and that the mean Galactic latitude for sources is significantly below the midplane: 
(^) = (-0.095 ±0.001)°. 

Subject headings: Galaxies:Milky Way — ISM: star formation 
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1. INTRODUCTION 

The advent of large-number bolometer arrays at good ob- 
serving sites has enabled the mapping of large sections of 
the sky in the millimeter continuum. These observations are 
ideally suited for unbiased searches for dense gas traced in 
dust thermal continuum emission. Even at the cold temper- 
atures and high densities found in molecular clouds, dust 
emission remains an optically thin tracer of the gas col- 
umn density. Despite some uncertainties in the properties 
of dust emission in these regimes, these observations rep- 
resent the most efficient way to detect high column den- 
sity features in the molecular interstellar medium (ISM). The 
original surveys of (sub)millimeter conti nuum emission fo- 
cused on n earby molecular clouds (e.g.. iMotte et al.l 119981 : 
iJohnstonFe t al. 2000; Enoch et al. 2007). With improvements 
in instrumentation and methodology, surveys of large sec- 
tions of the Galactic plane have become possible. This paper 
presents results from one such survey: the Bolocam Galactic 
Plane Survey (BGPS). The Bolocam instrument at the Cal- 
tech Submillimeter Observatory has been used to conduct 
a survey of the Galactic plane from -10° < £ < 90.5° and 
-0.5° < b < 0.5° in the continuum at A = LI mm. The sur- 
vey has been augmented with larger latitude coverage from 
77° < £ < 90.5° and studies of individual star forming re- 
gions near ^ = 1 1 1° and in the W3/4/5, Gem OBI and IC 1396 
regions. The observations, reduction and calibration of sur- 
vey da ta are presented in a companion paper by Aguir reet alJ 
(l2009h . hereafter Paper I. This paper presents the catalog of 
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millimeter continuum sources in the surveyed region. 

The interpretation of millimeter continuum maps from 
nearby molecular clouds is guided by preceding studies in an 
array of molecular gas tracers that developed a detailed under- 
standing of the density structure of clouds (see the review by 
lEvansll999l) . The millimeter continuum structures detected in 
the observations of nearby molecular clouds were molecular 
cores, i.e., dense structures that are linked to the formation of 
individual stellar systems. The detection of cores was driven 
by sensitivity to bright objects with angular scales of tens of 
arcseconds and filtering out larger angular scale structures. 

Although the sensitivity and filtering properties of the in- 
struments remain the same when directed at the Galactic 
plane, the typical distance to the objects becomes significantly 
larger (^^5 kpc vs. ^250 pc). In addition, the range of pos- 
sible distances spans foreground objects in nearby star form- 
ing regions to bright features viewed from across the Galac- 
tic disk. Hence, the physical nature of the objects seen in 
millimeter surveys of the Galactic plane is less clear than it 
is for nearby molecular clouds. Guided again by line stud- 
ies of nearby molecular clouds, surveys such as the BGPS 
are likely detecting a range of structures within the molecu- 
lar ISM. The actual objects that are detected are a product of 
matching the average column density of dust to the angular 
scales over which the survey is sensitive. Given the distance 
to the objects, many different types of objects may have an av- 
erage column density detectable on the ^ 30" scales to which 
the survey is most sensitive. 

In addition to the BGPS, other surveys are beginning to 
explore the Galactic plane in the millimeter dust contin- 
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uum. The ATLASGAL survey uses the Laboca instrument 
on the APEX telesco pe to survey the southern Galactic plane 
(lSchulleretal.ll2009h at A = 0.87 mm. After commission- 
ing, the SCUBA2 instrument will be used to survey the plane 
at A = 0.85 mm in the JCMT Galactic Plane Survey (JPS) 
(iMoore et al.l 120051) . These data will complement the open 
time Hi-GAL project of the Herschel Space Telescope which 
will ob serve the inner Galactic plane in several far infrared 
bands jNoriega-Crespo & Molinaril 120081) . The BGPS and 
these other surveys are conducted in a blind fashion, mean- 
ing that the survey spans all of a large region with no pref- 
erential targeting of subregions within that section. Since 
the survey regions span the low-latitude Galactic plane, they 
likely cover a large fraction of the star forming regions in the 
Galaxy. Hence, the maps will serve as the finding charts for 
studying dense gas in our Galaxy using the next generation of 
instrumentation, in particular the Atacama Large Millimeter 
Array. 

This paper presents a catalog of the BGPS data suitable 
for statistical analysis and observational follow-up of the mil- 
limeter objects in the Galactic plane. The catalog is produced 
by a tailor-made algorithm, Bolocat, which was designed and 
tested in concert with the data reduction and calibration of the 
BGPS images. Consequently, the algorithm specifically re- 
flects and accounts for the peculiarities in the BGPS images 
in a fashion that a generic source extraction algorithm could 
not. We present the catalog through a brief description of the 
BGPS data and Bolocat (® . We present our validation 
of the source extraction algorithm in ^and the catalog itself 
in ^ Using observations of the NH3 inversion transitions to- 
wards several of the objects, we are able to present constraints 
on the physical properties of catalog sources which are dis- 
cussed in ^ Finally, we present several summary products 
from the catalog in ^ 

2. THE BOLOCAM GALACTIC PLANE SURVEY 

The BGPS is described in full detail in Paper I. The sur- 
vey uses the 144-element Bolocam array on the Caltech Sub- 
millimeter Observatory, which observes in a band centered 
at 268 GHz (LI mm) and a width of 46 GHz (Glenn et al. 
[2003). Of particular note, the bandpass is designed to reject 
emission from the CO(2 — > 1) transition, which is the domi- 
nant line contributor at these wavelengths. The survey cov- 
ers 150 square degrees of the Galactic plane, primarily in the 
first Galactic quadrant. In the first quadrant, the survey ranges 
from -0.5° to 0.5° in latitude. Survey data are observed in per- 
pendicular scans made in the longitude and latitude direction, 
with the exception of one block of scans in the W5 region. All 
of the survey fields are observed in either l°xl°or3°xl° 
sub-fields which are then jointly reduced producing several 
individual maps of the plane spanning several degrees. 

Imaging of the time stream data requires correcting the sig- 
nal for the dominant atmospheric contribution in the bolome- 
ter array. Since the fundamental assumption in the correction 
is that the atmosphere is constant across the array, the data 
are not sensitive to spatial structures with sizes larger than 
the array (7.5'). More detailed analysis shows that structures 
larger than 3.5' are not fully recovered. An astrophysical im- 
age is generated by analyzing the variations of the signal in 
the time domain and identifying structures due to astrophysi- 
cal signal in the image domain. By iteratively refining the es- 
timates of instrumental and atmospheric contributions to the 
time stream, these effects can be decoupled from astrophys- 
ical signal producing a map containing only signal and irre- 



ducible noise. Of note for the cataloging process, the mapping 
procedure yields (1) the BGPS map (2) a model of the signal 
at each position or equivalently the residual of subtracting the 
model from the survey map and (3) a record of the number of 
times a bolometer sampled a given position on the sky. The 
pixel scale of the BGPS is set to 7.2" or 4.58 pixels across the 
33" FWHM beam. 

3. CATALOG METHODOLOGY 

The BGPS catalog algorithm was developed jointly with 
the mapping process in response to the goals of the BGPS. 
The BGPS is a large-scale survey of the Galactic plane in a 
comparatively unexplored waveband. The resulting images 
contain a wealth of structure as presented in Paper I. A sam- 
ple image appears in Figure[T]displaying many of the features 
that motivate the structure of the catalog algorithm. The map 
contains several bright, compact sources which frequently ap- 
pear in complexes. The bright sources commonly have low 
intensity skirts or filaments of fainter emission surrounding or 
connecting them. In addition, both bright and faint sources 
have resolved, irregular shapes. The brightest sources com- 
monly have negative bowls around them as a consequence of 
the PCA cleaning during the iterative mapping process (Paper 
I). These bowls can affect flux density measurements for these 
objects and depress the significance of surrounding emission. 
The maps also show noise with a spatially varying rms on 
scales from pixel (sub beam) to degree scales. Some artifacts 
that have not been completely removed by the imaging pro- 
cess still remain, manifesting as correlated noise in the scan 
directions. 

Given these features in the images, we opted to design 
our own cataloging algorithm rather than adapt a pree xisting 
routine to the BGPS data su ch as Clumpfind (William s et alj 
11994 or Source Extractor ("B ertin & Arnoutslll99 6). Since 
the physical nature of the objects is not well-constrained, we 
did not impose prior information on the objects that were to 
be extracted, such as the functional form of the intensity pro- 
file. Instead, the algorithm separates emission into objects 
defined by a single compact object and nearby emission. The 
catalog is driven by the need to recover and characterize the 
properties of the compact sources. The data retained in the 
catalog needed to be well-suited for follow-up observations, 
directing observers to the most interesting and bright features 
in the images. However, filamentary structure commonly ap- 
pears without an associated compact source and we devel- 
oped a routine that could also characterize isolated filaments 
of emission as objects. The routine also needed to account for 
the non-uniform nature of the noise properties over the maps. 
The algorithm described below is the result of applying these 
design criteria. 

We investigated the behavior of both Clumpfind and Source 
Extractor on a subset of the data to verify that a new approach 
was necessary. Our algorithm shares a common design philos- 
ophy with these other approaches. For the recommended user 
parameters, the Clumpfind algorithm tended to break regions 
into sub-beam sized objects based on small scale features in 
the objects from both real source structure and the spatially 
varying noise. After exploring the input parameters to Source 
Extractor, we were unable to satisfactorily identify filamen- 
tary structure in the images. Furthermore, those bright objects 
that were detected were not divided up into subcomponents. 
These shortcomings largely arose from (mis)applying these 
algorithms outside of the domain in which they were orig- 
inally developed. The unsatisfactory results suggested that 
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developing a tailored cataloging algorithm was necessary. 

Sources are extracted from the Bolocam images using a 
seeded watershed algorithm ( Yooii2004t) . This algorithm as- 
signs each pixel that is likely to contain emission to a single 
object based on the structure of the emission and nearby lo- 
cal maxima. Watershed algorithms are so named since they 
can be used to identify the watershed associated with a given 
body of water To make the analogy to the problem at hand, 
we consider the inverted emission map of a region as a two- 
dimensional surface. Several initial markers are identified (the 
seeds) associated with the lowest points of the inverted image 
(peaks in the original image). The surface is then "flooded" 
from the starting markers at successively higher levels. At 
each new level, newly flooded regions are linked to the origi- 
nal seeds associating each point on the surface with one of the 
starting seeds. 

Source identification takes place in a two-step process. 
First, significant regions of emission are identified based on 
the amplitu de o f the signal compared to a local estimate of 
the noise 03.1b . These sources are expanded to include ad- 
jacent, low-significance signal 03.2b . Second, each region so 
identified is examined and subdivided further into substruc- 
tures if ther e is s ignificant contrast between the sub-objects in 
the region ( ^3.31 ). 

The algorithm which accomplishes this object identification 
contains user-determined parameters. In the absence of phys- 
ical priors, we have set these parameters to reproduce a "by- 
eye" decomposition of some of the images. We use the same 
set of parameters to decompose all images in the BGPS and 
we have validated the parameter selection through in depth 
testing of the accuracy of the resulting catalog (®. Even 
though the choice of some parameter values are motivated by 
the desired behavior of the algorithm, the influence of these 
choices is well-known. 

3.1. Local Noise Estimation 

Objects in the BGPS catalog are identified based on their 
significance relative to the noise in the map. Owing to varia- 
tions in the observing conditions and integration time across 
an individual BGPS field, the underlying noise at each po- 
sition can vary substantially across a mosaic of several scan 
blocks. The spatially varying noise RMS (cr) is calculated as 
a function of position using the residual maps from the im- 
age production pipeline (Paper I). These maps consist of the 
actual images less the source emission model (where the sub- 
traction is performed in the time-stream). The residual maps 
are used in the map reconstruction process and are nominally 
source-free representations of the noise in the map. We esti- 
mate the RMS by calculating the median absolute deviation in 
an 11x11 pixel box (80" square) around each position in the 
map. The size of the box is chosen as a compromise between 
needing an adequate sample for the estimator and the need to 
reproduce sharp features in the noise map. The median abso- 
lute deviation estimates the RMS while being less sensitive to 
outliers relative to the standard deviation. The resulting map 
of the RMS shows significant pixel-to-pixel variation owing to 
uncertainty in the RMS estimator for finite numbers of data. 
To reduce these fluctuations, we smooth the RMS maps with a 
two-d imensional Savitsky-Golay smoothing filter (Press et al. 
Il992h . The filter is seventh degree and smoothes over a box 
12' on a side. The action of the filter is to reproduce large scale 
features in the map such as the boundaries between two noisy 
regions while smoothing noise regions on beam-size scales. 
An example map of the RMS is shown in Figure |2] after all 



processing (RMS estimate and smoothing). 

The statistical noise properties of the map are shown in Fig- 
ure [3] showing the intensity at each position in the map, / 
and the signal-to-noise ratio I /a. In general, generating an 
error estimate from the residual map tracks the position-to- 
position variation of the noise quite well. Transforming the 
data into significance space (i.e., dividing the data image by 
the noise estimate) yields a normal distribution with some de- 
viations. As seen in the bottom panel of the Figure in partic- 
ular, there is a wing of high significance outliers associated 
with the signal in the map. There are more negative outliers 
than expected, and this arises from pixel-to-pixel variations in 
the RMS based on local conditions. Since our method can- 
not reproduce these variations perfectly, regions with higher 
true RMS than represented in the RMS map will manifest as 
high significance peaks. There are likely regions for which 
the RMS is an overestimate and such regions contribute a 
slight excess of data near I/a = 0.0, which is impossible to 
discern in the face of the large number of data that are ex- 
pected there. We note the excess of negative outliers may also 
result from incomplete removal of astrophysical signal in the 
model map (see Paper I). In principle, it is possible to track 
the variations of the noise distribution on sub-beam scales us- 
ing information in the time domain. We found, however, that 
the generation of noise statistics from the time stream is unre- 
liable and produces non-normal distributions of significance. 
However, we do use the time-stream estimate of the noise to 
identify regions where the residual images are contaminated 
with astrophysical signal and we replace the RMS estimate 
at positions with spuriously large estimates of the noise with 
a local average. A Gaussian fit to the distribution of signif- 
icance values reproduces the shape of the distribution well 
with the exceptions noted above. Fits to the normal distribu- 
tion find the central value of the distribution (5') zero in all 
cases (-0.02 < {I /a) < 0.02). The width of the distribution, 
which should be identically 1 if a(£,b) is a perfect estimate 
of the noise, is usually within 2% of 1.00. Hence, the signif- 
icance distribution has been transformed to the best approxi- 
mation of a normal distribution with zero mean and dispersion 
equal to 1 . 

3.2. Identification of Significant Emission 

Given an estimate of the significance of the image at every 
position in the data, emission is identified as positive outlying 
data that are unlikely to be generated by noise. These data 
form connected regions of positive deviation above a given 
threshold. We initially mask the data to all regions above 
2(T(£,b). Owing to image artifacts from the Bolocam map- 
ping process, the emission profile shows low-amplitude noise 
on scales smaller than a beam. As a consequence of this noise, 
holes in the mask can appear where negative noise fluctua- 
tions drive a pixel value below the threshold. Additionally, 
a set cut at a given threshold tends to produce an artificially 
jagged source boundary. Finally, the PCA cleaning process 
creates negative bowls around strong sources. To reduce the 
influence of these effects, we perform a morphological open- 
ing operation on the data with a round stru cturing element 
with a diameter of a half a beam-width (see Doughertyl l 19921 
for details). The opening operation minimizes small-scale 
structure at the edges of the mask on scales less than that of a 
beam by smoothing the edges of the mask. 

Following the opening operation, the remaining regions are 
expanded to include all connected regions of emission above 
la(£,b) since regions of marginal significance adjacent to re- 
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Fig. 1 . — A sample subsection from the BGPS map near l = 24° where the units of the image are Jy beam ' . White contours indicate the shapes of features 
above the top end of the grayscale, beginning at 0.33 Jy beam"' and increasing by 0.33 Jy beam"' linearly. The field highlights several of the structures that 
informed the design of the cataloging algorithm, including bright compact sources and filamentary structures as well as spatially varying noise and negative bowls 
around strong sources. 



gions of emission are likely to be real. After this rejection, a 
morphological closing operation is performed with the same 
structuring element as the previous opening operation. The 
closing operation smoothes edges and eliminates holes in the 
mask. This processing minimizes the effects of the sharp 
cut made in significance on the mask, allowing neighboring, 
marginal-significance data to be included in the mask. An ex- 
ample of the effects for the opening and closing operations 
is shown in Figure |4] This results in our final mask contain- 
ing emission. To assess the effects of the mask processing on 
our data, we examined the brightness distribution of the pix- 
els in the mask pre- and post-processing by the opening and 
closing operations. Figure |5]presents the pixel brightness dis- 
tributions of the regions included in our mask both before and 
after applying the operators. The eliminated pixels are at low 
significance and the added pixels include marginal emission 
as desired. The processed mask also includes some negative 
intensity pixels in the objects. The added negative pixels are 
only 5% by number of the pixels added by mask processing 
a nd a re not included in the calculation of source properties 
(El. 

3.3. Decomposing Regions of Emission 



The emission mask is a binary array with the value 1 where 
the map likely contains emission and otherwise. Identifica- 
tion of objects occurs only within the emission mask region. 
We consider each spatially contiguous region of emission sep- 
arately and identify substructure within the region. We deter- 
mine whether this is significant substructure within the map 
by searching for distinct local maxima within the region. If 
there is only one distinct local maximum, the entire region is 
assigned to a single object. Otherwise, the emission is divided 
among the local maxima using the algorithm described below. 

To assess whether a region has multiple distinct local max- 
ima, we search for all local maxima in the region requiring 
a local maxima to be larger than all pixels within 1 beam 
FWHM of the candidate. For each pair of local maxima, we 
then identify the largest contour value that contains both local 
maxima ijcni) and measure the area uniquely associated with 
each local maximum. If either area is smaller than 3 pixels (as 
it would be for a noise spike) then the local maximum with the 
lowest amplitude is removed from consideration (see Figure 
|6]l. Furthermore, if either peak is less than 0.5cr above /c„, then 
the maximum with the smaller value of I^/a is rejected. This 
filter dramatically reduces substructure in regions of emission 
that is due to noise fluctuations. We note that the filter is fairly 
conservative compared to a by-eye decomposition. As a con- 
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Fig. 2. — Local variations in the noise for the W5 region of the BGPS. The 
noise map is generated from the residual map created in the image genera- 
tion process. The map illustrates variations in the noise level due to variable 
coverage by the survey and changing observing conditions. The W5 field 
was selected for display based on its varied survey coverage resulting in the 
observable structure in the noise map. 



sequence, compound objects are de-blended less aggressively 
than a manual decomposition. The rejection criteria are eval- 
uated for every pair of local maxima in a connected region, so 
every local maximum that remains after decimation is distinct 
from every other remaining maximum. 

The final set of local maxima are thus well-separated from 
their neighbors, have significant contrast with respect to the 
larger complex, and have a significant portion of the image 
associated with them. Using this set of local maxima, we par- 
tition the regions of emission into catalog objects. This parti- 
tioning is accomplished using a seeded watershed algorithm. 
Such an algori thm is, in principle, si milar to the methods used 
in Clumpfind jWilliams et al.ll 19941 ) in molecular l ine astron- 
omy or Source Extractor (Bertin & Arnoutsll 19961) in optical 
astronomy. To effect the decomposition, the emission is con- 
toured with a large number of levels (500 in this case). For 
each level beginning with the highest, the algorithm clips the 
emission at that level and examines all the resulting isolated 
sub-regions. For any pixels not considered at a higher level, 
the algorithm assigns those pixels to the local maximum to 
which they are connected by the shortest path contained en- 
tirely within the subregion. If there is only one local maxi- 
mum within the sub-region, the assignment of the newly ex- 
amined pixels is trivial. The test contour levels are marched 
progressively lower until all the pixels in the region are as- 
signed to a local maximum. An example of the partitioning of 
a region is shown in Figure |2l 

3.4. Measuring Source Properties 

A BGPS catalog object is defined by the assignment al- 
gorithm established above which yields a set of pixels with 
known Galactic coordinates i£j,bj) and associated intensities 
Ij (Jy beam"'). Source properties are determined by emission- 
weighted moments over the coordinate axes of the BGPS im- 
age. Positions in the mask with negative intensity are not in- 
cluded in the calculation of properties. For example, the cen- 



troid Galactic longitude of an object is 

with been defined in an analogous fashion. The sizes of the 
objects are determined from the second moments of the emis- 
sion along the coordinate axes: 



T.i{bj-bce,?j^Ij 



(2) 
(3) 
(4) 



The principal axes of the flux density distribution are deter- 
mined by diagonalizing the tensor 



1 = 



(5) 



The diagonalization yields the major {a„aj) and minor ia„,i„) 
axis dispersions as well as the position angle of the source. 
Since the moments are calculated for and emission mask that 
has been clipped at a positive significance level, these esti- 
mates will slightly underestimate the sizes of the sources lead- 
ing to size estimates that could be smaller than the beam size. 
Once projected into the principal axes of the intensity distri- 
bution, the angular radius of the object can be calculated as the 
geometric mean of the deconvolved major and minor axes: 



Or = V [(cr; 



2 

ma J 



,1/4 



(6) 



Here abm is the rms size of the beam (= 6'FWHM/\/81n2 = 14") 
and 77 is a factor that relates the rms size of the emission dis- 
tribution to the angular radius of the object determined. The 
appropriate value of 77 to be used depends on the true emis- 
sion distribution of the object and its size relative to the beam. 
We have elected to use a value of 77 = 2.4 which is the median 
value derived for a range of models consisting of a spherical, 
emissivity distribution ae(r) cx r'* with -2 < 7 < 0, a range 
of radii relative to the beam 0.1 < Or/Oy^hu < 3 at a range 
of significance (peak signal-to-noise ranging from 3 to 100). 
The range of -q for the simulations is large, varying by more 
than a factor of two. If a specific source model is required for 
a follow-up study, the appropriate value of -q should be calcu- 
lated for this model and the catalog values can be rescaled to 
match. 

Uncertainties in the derived quantities are calculated by 
propagating the pixel-wise estimates of the flux density un- 
certainties through the definitions of the quantities. The co- 
ordinate axes are assumed to have negligible errors. The mo- 
ment descriptions of position, size and orientation are exact 
for Gaussian emission profiles and generate good Gaussian 
approximations for more complex profiles. An example of 
the derived shapes of the catalog objects is shown in Figure]?] 
for comparison with the shapes of the regions from which they 
are derived. In comparing the images, it is clear that the irreg- 
ularly shaped catalog objects can produce elliptical sources 
that vary from what a Gaussian fit would produce. However, 
in most cases the moment approximations characterize the ob- 
ject shapes well. 
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Fig. 3 . — Noise properties for the W5 BGPS field. The top panel shows the histogram of raw data values present in the image, in particular showing a large tail 
towards high values which arise because of both signal and high noise values at the map edge. The middle panel shows the significance distribution of the image 
after normalization by the noise estimate derived from the weight map. The bottom panel repeats the data from the middle panel but using a logarithmic scale 
for the y-axis, highlighting the behavior of the distribution at large values of the significance. In both the middle and the bottom panel, a gray curve shows the 
Gaussian fit to the significance distribution. 



Since the extracted objects are often asymmetric, the cen- 
troid of the emission does not correspond to the peak of emis- 
sion in an object. Because of this possible offset, we also 
report the coordinates of the maximum of emission. To min- 
imize the effect of noise on a sub-beam scale, we replace the 
map values with the median value in a 5 x 5 pixel (36") box 
around the position and determine the position of the maxi- 
mum value in the smoothed map. While this reduces the po- 
sitional accuracy of the derived maximum, it dramatically im- 
proves the stability of the estimate. Since the peaks of emis- 
sion are rarely point sources, the loss of resolution is made up 
for by the increased signal-to-noise in the map. These coor- 
dinates are reported as £„,ai and bi„ax- If the catalog is being 
used to follow-up on bright millimeter sources in the plane 
these coordinates are preferred over the centroid coordinates. 
Moreover, since the maxima are less subject to decomposi- 



tion ambiguities, we name our sources based on these coor- 
dinates: GMilM ± bb.bbb. The smoothed map is only used 
for the calculation of the position for the maximum; all other 
properties are calculated from the original BGPS map. 

The flux densities of the extracted objects are determined 
by measuring the flux density in an aperture and by inte- 
grating the flux density associated with the region defined 
by each object in the watershed extraction. In the catalog, 
three flux density values for each source are reported, corre- 
sponding to summing the intensity in apertures with diameters 
dapp = 40", 80", and 120". We refer to fliese flux densities 
as 540,580, and 5i2o respectively. Because many of the ob- 
jects are identified in crowded fields, we refrain from a gen- 
eral "sky" correction since inspection of the noise statistics 
suggests that the average value of the sky is quite close to 
zero. One notable exception to this case is where there are 
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Original Mask 




Mask After Closure 



Fig. 4. — Results of applying a morphological opening and closing operations to a mask. The left-hand panel shows the results of a clip at 2a and the right-hand 
panel shows the mask after it has been processed by the morphological operations. The opening eliminates artificial structure on scales smaller than that of the 
beam. The closure smoothes edges and eliminates holes. This largely cosmetic process produces smooth features which enables better source decomposition. 
The sample region is 24' on a side. 
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Added in Processing 
Removed in Processing 




0.2 0.4 O.f 

Intensity (Jy beam" 



Fig. 5. — The brightness distribution of pixels included in the emission 
mask pre- and post- processing by the opening and closing operators. The 
added pixels are mostly positive low significance pixels (dashed line) and the 
removed pixels (gray line) are confined to relatively faint emission. A small 
fraction (5%) of the added pixels have negative intensities. 



pronounced negative bowls in the emission around the object. 
Because large apertures may contain parts of these negative 
bowls, the aperture flux densities of bright objects are likely 
underestimated. In addition, we also report the object flux 
density value iS), which is determined by integrating the flux 
density from all pixels contributing to the object in the water- 
shed decomposition of the emission: 



(7) 



where A57, is the angular size of a pixel. The object flux 
density cannot, by definition of the regions, contain negative 
bowls, but the bowls indicate that emission is not fully re- 
stored to these values, which are thus likely underestimates 
of the true flux density for bright objects. Since objects typi- 
cally have angular radii ^ 60", 120" apertures do not contain 
all the flux density. Hence, S > 5i2o > ^so > ^40. Since the 




Fig. 6. — Schematic contour diagram of elimination criterion for local max- 
ima. Each of a pair of local maximum (here I\ , I2) has an area uniquely as- 
sociated with that local maximum (i.e., above a critical contour /«»). If the 
difference between /, /ct - /cr is smaller than 0.5 then the smaller ampli- 
tude local maximum is removed from the list of candidate local maxima. If 
the smaller area uniquely associated with a maximum (here A2) is smaller 
than 3 pixels, then that maximum is removed. 



surface brightness at e ach position in the map has an associ- 
ated uncertainty ( ^3. It , we calculate the the uncertainty in the 
flux density measurements by summing the uncertainties in 
quadrature. We further account for the uncertainty in the size 
of the beam owing to the uncertainties of the bolometer po- 
sitions in the focal plane of the instrument (see Paper I). The 
latter effect produces a fractional uncertainty in the beam area 
of d^lbm/^bm = 0.06. Thus, the uncertainty in a flux density 
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Fig. 7. — Example of the watershed algorithm apphed to the same region shown in Figure|4]('e/i') and the measured source shapes (right). The grayscale image 
shows the actual BGPS map image which runs on a linear scale from to 0.2 Jy beam"' . In the left panel, the blue boundaries indicate the regions identified 
using the seeded watershed algorithm. In the right panel, the ellipses indicate the shapes derived from the moment methods applied to the extracted objects. The 
example region shown is from the oiiginal £ = 32° BGPS maps and is 24' on a side. 



measurement (either aperture or integrated) is: 




where the ratio 57,/51i„„ accounts for the muhiple pixels per 
independent resolution element. The flux density errors do 
not account for any systematic uncertainty in the calibration of 
the image data which is estimated to be less than 10% (P aper 
I). We discuss the flux densities in more detail in Section [7721 
below. 

4. SOURCE RECOVERY EXPERIMENTS 

In order to assess the properties of catalog sources com- 
pared to the actual distribution of emission on the sky, we have 
created simulated maps of the emission, inserting sources 
with known properties. The generation of these simulated im- 
ages relies on the actual observed data for the £ = 111° re- 
gion of the survey. This field was selected for analysis based 
on its uniform noise properties and well-recovered astrophys- 
ical signal. Most of the astrophysical signal is identified in 
the time series data from the instrument and removed, leav- 
ing only irreducible (i.e., shot) noise in the time series with 
small contributions due to atmospheric, instrumental (e.g., the 
pickup of vibrations) and astrophysical signals (see Paper I for 
more detail). Then a series of Gaussian objects was sampled 
into the data time stream. The simulation time streams are 
then processed in a fashion identical to the final BGPS survey 
images resulting in a suite of test images to test the catalog 
algorithm on. Gaussian brightness profiles do not adequately 
represent all the structure seen in the BGPS maps, but they do 
capture the essential features of compact sources. The com- 
plicated interplay of extended structure with the mapping pro- 
cess likely affects the validity of the following tests. However, 
the behavior for point-like or small sources should be well re- 
produced. Finally, some of the astrophysical signal remains in 
the residual image which affects some of our tests to a small 
degree. 



4.1. Completeness 

The primary test of the algorithm is the degree of complete- 
ness of the catalog at various flux density limits. To this end, 
we have inserted a set of point sources into the time stream 
with flux densities uniformly distributed between 0.01 Jy to 
0.2 Jy arranged in a grid. We then process the simulated data 
with the source extraction algorithm and compare the input 
and output source catalog. The results of the completeness 
study are shown in Figure [8] For this field, the mean noise 
level is (T = 24 mJy beam"'. The catalog is > 99% complete 
at the 5(7 level. Below this level, the completeness fraction 
rolls off until no sources are detected at 2(7. For this regular 
distribution of sources, there are no false detections through- 
out the catalog (even at the 2a level). The catalog parame- 
ters described above are selected to minimize the number of 
false detections. Even though the thresholding of the map is 
performed at the relatively low 2.0a level, the wealth of re- 
strictions placed on the extracted regions combine to raise the 
effective completeness limit to 5a. We note that the grid ar- 
rangement of test sources means that confusion effects are not 
accounted for in the flux limit. 

In Figure|9l we present the estimated completeness limit as 
a function of Galactic longitude for the BGPS. The ^=111° 
region has a relatively low value for the noise RMS with sev- 
eral fields being up to a factor of 4 larger than the complete- 
ness test field. We assume that the completeness limit scales 
with the local RMS. The assumption is motivated since all cat- 
alog source extraction and decomposition is performed in the 
significance map rather than in the image units. Thus, the ab- 
solute value of the image data should not matter Over > 98% 
of the cataloged area, the completeness limit is everywhere 
less than 0.4 Jy, thus the survey, as a whole can be taken to be 
98% complete at the 0.4 Jy level. 

4.2. Decomposition Fidelity 

In addition to completeness limits, we investigated the fi- 
delity of the algorithm by comparing extracted source proper- 
ties to the input source distributions. We compared the ex- 
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Fig. 8. — Completeness fraction deiived from fake source tests in a survey 
field. The fraction of sources recovered is plotted as a function of input source 
flux density. The vertical dotted lines indicate { 1,2, 3,4,5 )cr. The complete- 
ness tests indicate that the catalog is complete at the > 99% limit for source 
flux densities > 5cr. 

traded flux density values to the input flux density values 
as shown in Figure [TO] These simulations compare input 
and recovered properties for resolved sources with intrinsic 
FWHMs equal to 2.4 times that of the beam. The flux den- 
sities of the sources are varied over a significant range, and 
we measure the recovered values for the flux density using 
the 40", 80" and 120" apertures and we also compute the in- 
tegrated flux density in the object. Since the small apertures 
only contain a fraction of the total flux density, their average 
flux density recovery is naturally smaller than the input. As 
the aperture becomes better matched to the source size, the 
flux density recovery improves, in agreement with analytic 
expectations based on the source profile. Since many objects 
are larger than the 40" and 80" apertures, catalog flux density 
values can be viewed in terms of surface brightnesses within 
these apertures. 

We completed an additional suite of simulations, including 
sources with sizes significantly larger than the beam. These 
simulations are used in Paper I to test the effect of the de- 
convolution process, but we also examine the results of the 
source extraction. In Figure[TO]we show the recovered source 
size as a function of input size. In general, smaller sources are 
better recovered with objects much larger than the beam hav- 
ing a smaller size found. Radii of objects become biased for 
sources with sizes > 200", consistent with the experiments 
on deconvolution discussed in Paper I. Because of the im- 
age processing methods used, objects with sizes larger than 
this should be treated with caution. However, we note that no 
sources in the actual catalog are found with sizes > 220" and 
the mean size is only 60" so objects of these sizes and their 
properties are likely to be recovered. The catalog algorithm 
thus does not introduce any additional bias into the proper- 
ties of the recovered sources. However, the sizes of objects 
may not be completely recovered in the mapping process and 
thus the extracted properties may not accurately reflect the 
true properties of sources because of limitations in the map 
making process. These map properties should be viewed with 
the caveats presented in Paper I. 

4.3. Deblending 



A primary question of the catalog algorithm is how well 
it separates individual sources into their appropriate subcom- 
ponents. Any approach to source extraction and decomposi- 
tion will be unable to separate some source pairs accurately. 
However, the limiting separation for source extraction is de- 
termined by the method. Our approach is further complicated 
by the assignment of emission from every position in the cat- 
alog to at most one object. Thus, blended sources are diffi- 
cult to extract separately unless the individual local maxima 
t hat d efine the sources are apparent according to the criteria in 
^3.31 To assess how the catalog algorithm recovers sources, 
we generated maps of sources with a random distribution in 
position and processed the maps with the catalog algorithm. 
The input sources have elliptical Gaussian profiles with a uni- 
form distribution of major and minor axis sizes between our 
beam size (33") and 72", which represent typical objects seen 
in the BGPS maps. The input sources are oriented randomly. 
For each pair of input sources, we determine whether the pair 
is assigned to the same catalog object. We then measure the 
fraction of blended sources as a function of separation. The 
results of the experiment are shown in Figure [TT] which indi- 
cates that the algorithm resolves most sources that are sepa- 
rated by > 75" or 2 beam FWHM. For point sources, the nom- 
inal resolution limit is w 1 beam FWHM; however the algo- 
rithm requires ^2 FWHM between local maxima so this limit 
on deblending is consistent with the design. Furthermore, 
the test objects are resolved, which complicates the analysis. 
Thus, Figure[TT]also indicates the fraction of blended sources 
for pairs of the smallest sources and for pairs of the largest 
sources (i.e., both sources are in the bottom 15% of the size 
distribution or the top 15% respectively). In separating the 
data based on size, we can assess the degree to which source 
structure will complicate deblending. As expected, pairs of 
small sources are easier to resolve (50% are separated for dis- 
tances of 60") and pairs of large sources more difficult to 
deblend (50% at 85"). The large source curve crosses the 
population curve at large separations (120") because of small 
number statistics rather than any algorithmic feature. 

5. THE CATALOG 

The BGPS catalog contains 8358 sources extracted using 
the algorithm described in ^ An excerpt from the catalog is 
shown in Table [T] The full catalog and the BGPS image data 
are hosted at the Infrared Processing and Analysis Center'. 

6. THE PHYSICAL PROPERTIES OF CATALOG 
SOURCES 

The BGPS data contain a wealth of sources extracted from 
the images with a wide variety of source properties. The 
physical meaning of a catalog source remains unclear, how- 
ever. Ascribing nomenclature such as a "clump" or a "core" 
to catalog objects may be misleading in certain cases as these 
terms can refer to specific density, mass and size regimes (e.g., 
IWilliams et al.l l2b00). In this section, we explore the range of 
parameter space a BGPS catalog object may occupy. We also 
refer readers to two additional discussions of BGPS sources 
in the Galactic center (Bally et al. 2009) and in the Gem OB I 
star forming region in the outer Galaxy (iDunham et al.l l2009'). 
These two studies consider the source properties in regions 
where most of the sources are likely found at a known dis- 
tance. In contrast, blind observations of the Galactic plane 
can span a range of distances. 

' |http : //irsa ■ ipac ■ caltech . edu/data/BOLOCAM_GPS/| 
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Fig. 9. — The 99% completeness limit in flux density for point som'ces plotted as a function of Galactic longitude. The completeness limit is 5 times the median 
RMS value in 1 ° bins across the survey. The significant variation in the completeness limit stems from the variable observing conditions and integration time for 
each field in the survey. 
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Fig. 10. — Recovery of source flux density and source size for simulated observations. The left-hand panel compares input and recovered source flux densities 
for resolved objects (FWHM = lAOtgam = 79"). The right-hand panel shows the size measured in the simulated objects for an array of input object sizes. The 
vertical dashed Une at 210" indicates the largest major-axis size of objects actually found in the BGPS implying source sizes are typically well recovered. 



To facilitate understanding BGPS source properties, we in- 
clude Robert F. Byrd Green Bank Telescope (GBT) observa- 
tions of the NH3(1,1) and (2,2) inversion transitions toward a 
set of BGPS detections in the ^ ^ 3 1 ° field. In making this 
comparison, we assume that the ammonia emission comes 
from the same regions as are responsible for the millime- 
ter continuum emission, an assum ption supported by stud- 
ies in local molecu lar clouds (e.g.. iRosolow skv et al.l 120081 : 
iFriesen et alj 120091) . These data are presente d formally and 
in much greater depth elsewhere (iDunham et al.ll2009l) . Here 



we note that the ammonia data provide complementary infor- 
mation on the line-of-sight velocities (Vlsr), line widths ay 
and kinetic temperatures of dense gas Tk- The velocity infor- 
mation, in particular, yields distance estimates to the sources 
by assuming that the radial motion of the line emitters arises 
entirely from their orbital motion ar ound the Galactic center. 
We then use the kinematic model of iReid et alJ ( 120091) to es- 
timate the distances to BGPS objects. The Reid et al. model 
uses VLBl parallaxes to establish the distances to several high 
mass star forming regions. One of the primary results of their 
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TABLE 1 
The BGPS Catalog (Excerpt)^ 



No. 


Name 


^max 




e 


b 






PA 


Or 


540 


^80 


^120 


5 






(°) 


(°) 


(°) 


(°) 


(") 


(") 


(°) 


(") 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


1 


G000.000 + 00.057 


0.000 


0.057 


0.004 


0.063 


32 


17 


89 


38 


0.19±0.06 


0.56±0.11 


1.01±0.18 


0.67±0.12 


2 


G000.004 + 00.277 


0.004 


0.277 


0.007 


0.276 


39 


34 


145 


80 


0.21±0.04 


0.59±0.09 


0.99±0.13 


1.40±0.17 


3 


GOOO.006-00.135 


0.006 


-0.135 


0.019 


-0.138 


43 


32 


115 


82 


0.22±0.05 


0.71 ±0.11 


1.39±0.17 


1.41±0.19 


4 


GOOO.010 + 00.157 


0.010 


0.157 


0.019 


0.156 


62 


24 


80 


81 


0.61±0.06 


1.62±0.15 


2.33±0.21 


3.55±0.30 


5 


GOOO.016-00.017 


0.016 


-0.017 


0.012 


-0.014 


50 


35 


44 


93 


1.24±0.10 


3.68±0.26 


6.23±0.43 


10.94±0.74 


6 


GOOO.018-00.431 


0.018 


-0.431 


0.016 


-0.431 


23 


8 


97 




0.08±0.05 


0.16±0.09 


0.19±0.14 


0.15±0.07 


7 


G000.020 + 00.033 


0.020 


0.033 


0.019 


0.036 


38 


32 


54 


77 


0.83±0.09 


2.68±0.22 


4.86±0.37 


6.07±0.46 


8 


G000.020 - 00.051 


0.020 


-0.051 


0.021 


-0.052 


34 


26 


85 


62 


1.14±0.10 


3.43±0.25 


5.75±0.41 


6.37±0.46 


9 


G000.022 + 00.251 


0.022 


0.251 


0.025 


0.250 


26 


14 


90 




0.10±0.04 


0.24±0.09 


0.43±0.13 


0.26±0.08 


10 


G000.034- 00.437 


0.034 


-0.437 


0.037 


-0.438 


17 


11 


92 




0.06±0.05 


0.14±0.09 


0.18±0.14 


0.11 ±0.06 


11 


G000.042 + 00.205 


0.042 


0.205 


0.041 


0.208 


44 


24 


155 


70 


0.36±0.05 


1.21±0.12 


2.15±0.19 


2.61±0.23 


12 


GOOO.044-00.133 


0.044 


-0.133 


0.049 


-0.135 


29 


20 


97 


45 


0.14±0.05 


0.42±0.10 


0.76±0.16 


0.62±0.12 


13 


G000.044- 00.285 


0.044 


-0.285 


0.051 


-0.282 


26 


10 


137 




0.08±0.04 


0.33±0.08 


0.60±0.13 


0.15±0.06 


14 


GOOO.046 + 00.111 


0.046 


0. 1 1 1 


0.049 


0.108 


39 


24 


142 


64 


0.33±0.05 


0.92±0.11 


1.40±0.16 


1.83±0.20 


15 


GOOO.046-00.041 


0.046 


-0.041 


0.047 


-0.037 


34 


29 


168 


67 


0.47±0.07 


1.67±0.17 


3.36±0.28 


3.21±0.28 


16 


G000.046- 00.297 


0.046 


-0.297 


0.043 


-0.294 


30 


21 


141 


49 


0.11±0.04 


0.32±0.08 


0.52±0.12 


0.48±0.10 


17 


G000.048 + 00.065 


0.048 


0.065 


0.047 


0.070 


13 


12 


127 




0.14±0.05 


0.56±0.11 


1.34±0.17 


0.19±0.06 


18 


G000.052 + 00.027 


0.052 


0.027 


0.048 


0.029 


39 


33 


161 


79 


1.30±0.11 


3.78±0.27 


6.05±0.43 


8.10±0.56 


19 


GOOO.052-00.159 


0.052 


-0.159 


0.051 


-0.159 


21 


19 


48 


33 


0.11±0.05 


0.33±0.09 


0.67±0.14 


0.33±0.09 


20 


G000.054- 00.209 


0.054 


-0.209 


0.055 


-0.207 


37 


36 


73 


80 


0.87±0.07 


2.56±0.18 


4.09±0.29 


5.76±0.41 


21 


G000.066 + 00.209 


0.066 


0.209 


0.067 


0.211 


34 


28 


115 


66 


0.31±0.04 


l.OOiO.lO 


1.80±0.16 


1.86±0.17 


22 


G000.066 - 00.079 


0.066 


-0.079 


0.068 


-0.076 


41 


31 


106 


79 


2.20±0.16 


6.92±0.45 


11.24±0.72 


15.98±1.03 


23 


GOOO.068 + 00.241 


0.068 


0.241 


0.067 


0.244 


26 


17 


52 


37 


0.09±0.04 


0.25±0.08 


0.37±0.12 


0.34±0.09 


24 


GOOO.070 + 00.175 


0.070 


0.175 


0.075 


0.173 


47 


38 


159 


94 


0.49±0.05 


1.52±0.12 


2.71±0.20 


4.86±0.34 


25 


G000.070- 00.037 


0.070 


-0.037 


0.071 


-0.035 


27 


20 


36 


43 


0.29±0.06 


0.85±0.13 


1.58±0.21 


1.17±0.17 



Note. — (1) Running Source Number. (2) Name derived from Galactic coordinates of tiie maximum intensity in the object. (3)-(4) Galactic coordinates of maximum intensity 
in the catalog object. (5)-(6) Galactic coordinates of emission centroid. (7)-(9) Major and minor axis 1/f widths and position angle of source. (10) Deconvolved angular size of 
source. (1 1)-(13) Flux densities derived for 40", 80" and 120" apertures. (14) Integrated flux density in the object. The derivation of these quantities is described in ^3.41 

sources and ammonia observations. We match each ammonia 
spectrum to the BGPS catalog source for which the pointing 
center for the GET observation falls into the region assigned 
to the BGPS source. If multiple GET observations fall within 
a EGPS source, we consider only the spectrum with the high- 
est integrated intensity in the (1,1) transition. In this fashion, 
we match 42 ammonia spectra with EGPS sources. We esti- 
mate the line properties of the ammonia spectra (i.e., Vlsr, c^v, 
Tk) using the model of Rosolowskv et al. (2008). Of partic- 
ular note, there is only one velocity component of ammonia 
emission along each line of sight examined, supporting our 
association between EGPS emission and ammonia emission 
for these sources. 

Using the LSR velocities of the EGPS sources, we compare 
these objects to the broader distribution of molecular emission 
in this field. In particular, we compare where these sources are 
found to the "CO (1-^0) data of iJackson et"!!] ( 120061) . We 
present a latitude-integrated map of the molecular emission 
for the 1=3)1° region in Figure [12] Of note in this figure, all 
EGPS sources are associated with the '^CO emission at the 
> 2(7 level. However, the EGPS sources are not necessarily 
found at the brightest positions in the '^CO data cube. 

Towards these Galactic latitudes, the '^CO emission is 
found predominantly near Vlsr = 40 km s"' and 90 km s"'. 
EGPS sources are associated with both velocity features sug- 
gesting a wide range of kinematic distances. We explore the 
association of EGPS emission with the tracers of low den- 
sity molecular gas in the full presentation of these results 
(Dunham et al. 2009). 

To explore the properties of EGPS sources in more detail, 
we estimate the source masses from the flux density of the 
continuum emission combined with distance and temperature 
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Fig. 1 1. — Fraction of source pairs that remain blended as a function of 
source pair separation. The solid line shows the blending fraction for all 
source pairs. Since the sources are elhptical Gaussians, the dashed and dotted 
curves show the same function plotted only for the source pairs with mean 
sizes in the top and bottom 15%, indicating that large resolved sources are 
more difficult to separate algorithmically. 



work is a differentiation between the kinematic properties of 
high mass star forming regions and other components of the 
Galactic disk. Since the EGPS sources are likely associated 
with (high-mass) star forming regions, their refined model is 
appropriate for our work. 

The ammonia observations were taken before the EGPS 
catalog was complete, and the GET targets were selected by- 
eye from an early version of the EGPS images. Hence, there 
is not a one-to-one correspondence between EGPS catalog 
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Fig . 12. — Latitude-integrated '""CO ( 1 0) emission with the (l, V) loca- 
tions of 42 BGPS/NH3 sources indicated. All BGPS sources are co-located 
in the data cube with significant '^CO emission, though not necessarily the 
brightest locations. 



estimates from the GBT data: 



M-- 



:13.1M0 



1 kpc 



IJy 



exp(13.0 K/r)-l 
exp(l 3.0/20)-! 



(9) 



10) 



Here, we use the total flux density in the object S measured 
in the catalog. This estimate uses a T = 20 K dust temper- 
ature and averages the emergent spectrum across the Bolo- 
cam band assuming the dust emissivity scales like A~'^ and an 
opacity of k = 0.01 14 cm^ g"' (lEnocti et al.ll2006h . For these 
assumptions, the mean frequency of the Bolocam passband is 
v = 21\.\ GHz. We plot the masses of the BGPS sources iden- 
tified as a function of line-of-sight distance in Figure [13] For 
source temperatures, we use the kinetic gas temperatures de- 
rived from the ammonia observations. We use kinematic dis- 
tances from the LSR velocity in the ammonia observations. 
Towards these Galactic longitudes, the velocity-distance re- 
lation is double valued resulting in the distance ambiguity. 
We plot the mass and distances derived for each source at 
both possible distances; so each point appears twice in the 
plot, once on each side of the d = 7.1 kpc line. This 7.1 kpc 
distance is determined by the tangent velocity at these lon- 
gitudes. We also plot a curve representing the derived mass 
of the 5(7 completeness limit for point sources in this field 
(where = 0.14 Jy beam"') with a temperature of T = 20 K. 
We note that sources are not distributed uniformly down to the 
completeness limit since the ammonia observing targets were 
selected by eye and preferentially selected bright objects. 

Figure [13] illustrates the wide range in possible masses ex- 
tracted for the sources extracted from the survey data. First, 



depending on the resolution of the kinematic distance ambi- 
guity, the masses in this field can vary by nearly an order 
of magnitude. Even restricting to near distances, the mass 
range of detected objects spans over two orders of magnitude. 
While a 15 Mq object at the low mass range may indeed be a 
high mass core, the larger mass scales at farther distances are 
more appropriately linked to clumps, thos e large structures 
which form stellar clusters (IWilliams et al.ll2000.) . This range 
is largely the result of the wide range of object flux densities 
recovered from the original maps (see the right-hand panel 
of Figure [T3Tl. but naturally also depends strongly on the true 
distance to the source. 

In Figure [14] we plot the derived values of the decon- 
volved projected radius R = d Or and particle density ne = 
3M / (4ttR^ fiiriii) as functions of distance. We have taken the 
mean particle mass to be /i = 2.37 (Kauffmannet al. 2008). 
In Figures [13] and [14] we plot the behavior of sources as a 
function of distance to illustrate the importance of distance 
determinations in describing the physical properties of these 
sources. In these plots, dotted lines indicate roughly the lim- 
iting sizes and densities to which the survey is sensitive. Be- 
cause the ability to estimate the deconvolved source sizes de- 
pends on the signal-to-noise of the observations, the size limit 
9r = 33"(= 1 FWHM) is chosen as a guide rather than a hard 
limit. It is formally possible resolve sources with sizes smaller 
than the beam size (Equation O though the stability of such 
deconvolutions is suspect for Or < Ot,,,. 

These plots emphasize the difference in objects recovered 
by the BGPS observations of the Galactic plane compared 
to observations of nearby molecular clouds where the Bolo- 
cam beam and sensitivit y is well-match e d to the size of a 
star-forming core (e.g. Youn g et al.ll2006l: lEnoch et al.ll2006l 
120071) . While the BGPS catalog contains sources with core- 
like properties (n > 10^^ cm"-', R <Q.l pc), most of the re- 
covered sources have larger sizes and lower density scales. 
The projected mass and size limits [Figure [T3] (left) and Fig- 
ure [14] (left)] limit the detection of individual star-forming 
cores in the BGPS to objects within d < I kpc. Beyond this 
distance, BGPS is limited to the detection of clumps within 
larger clouds. As is illustrated by Figure[T4](right), the density 
of gas to which the survey is sensitive, drops off precipitously 
with distance to the emitting source. Objects in the catalog 
are thus a range of clumps and cores, and their detection de- 
pends on their distance, physical size, and dust temperature. 
Despite these caveats, we note that most sources in this field 
have mass and density scales typical of massive, dense clumps 
in local clouds. Given the link that these clumps have to local 
cluster formation, it is likely that most sources in the BGPS 
are associated with the formation of stellar clusters. 

7. STATISTICAL PROPERTIES OF BGPS SOURCES 

In this section, we describe the basic statistical representa- 
tions of the source catalog. In all cases, the reader is cautioned 
that the source catalog is derived from the BGPS maps which 
are a filtered representation of the true sky brightness at any 
given position. In particular, large scale emission is filtered 
out of the final maps on scales larger than 3.5' (Paper I). 

7.1. The Galactic Distribution of Flux Density 

Figure [15] presents the total flux density extracted from the 
images in catalog sources as a function of coordinates in the 
Galactic plane. In total, the catalog sources contain 11.9 kJy 
of flux density. Difficulties in recovery of large scale struc- 
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Fig. 13. — Mass estimates for BGPS sources as a function of line-of-sight distance (left) and flux density (right). The masses are derived from a combination 
of the BGPS data (for dust continuum flux densities) and ammonia data (for gas temperatures and kinematic distances) for matched sources in the £ = 32° field. 
This figure highlights the wide range in the masses of objects recovered in the BGPS catalog. 
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Fig. 14. — Radius (left) and density (right) estimates for BGPS catalog objects matched with ammonia observations in the £ = 32° field. Like Figure [Tsl a 
given object is plotted twice for possible near and far kinematic distance estimates. The dotted line indicates the projected size for 9r = 33" (1 beam FWHM) in 
the left-hand panel. In the right-hand panel, the dotted line indicates the behavior for a source with S = 0.15 Jy =5cr and Br = 33". Sources can appear below the 
6r = 33" line because the beam deconvolution (Equation|6). 



ture notwithstanding, this represents the emission profile of 
the Galaxy along the longitude and latitude direction. 

The right-hand panel of Figure [15] shows that the distribu- 
tion of flux density in sources peaks at the flux-weighted mean 
(b) = -0.095 ±0.001 degrees. Such an offset was also seen in 
the ATLASGAL survey of t he 12° < ^? < -30° r egion of the 
Galactic plane at A=870 fim (S chuUer et al.l20 ()^. They spec- 
ulated that the offset may result from the Sun being located 
slightly above the Galactic plane. To further investigate this 
possibility, we examined the latitudinal flux density distribu- 
tion of BGPS sources as a function of longitude and compared 
these results to other tracers of the Galactic disk. In particu- 
lar, we compared the flux density weighted values of the mean 
latitude for objects in our c atalog to the COBE/DIRBE data 
at A = 2.2 /xm and 200 ^m dHauseret al.lll998) and the inte- 
grated '^CO (1-0) data of Dame et al. (2001). The results of 
the comparison are shown in Figure [16] For all data sets, we 
only compute the flux-density weighted average latitude over 



the range -0.5° < b < 0.5° which is the latitude range over 
which the BGPS is spatially complete. The data are averaged 
in 10° bins. The offset of BGPS sources appears at nearly 
all Galactic longitudes, but is most pronounced towards the 
Galactic center (-10° <£< 25°). The CO (1-0) data follow 
a similar trend with a significant offset towards negative lati- 
tudes. We verified that the effects persist in the CO data if the 
averaging range is expanded to larger latitudes (foreground 
emission from the Aquila rift affects the averaging for large 
latitude ranges \b\ < 5°). We also note that the emission pro- 
files from the DIRBE data do not follow the offset trends seen 
in the millimeter continuum or the CO. The DIRBE 2.2 fim 
data primarily traces emission from stellar photospheres and 
is affected by dust extinction and the 200 fim traces warm 
dust emission. Both of these emission features show a more 
symmetrical distribution around the Galactic midplane than 
do those tracers associated with the molecular gas. We con- 
clude that the offsets towards negative latitude are primarily 
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associated with the molecular ISM. The significant variations 
with the mean offset are likely the result of individual star 
forming structures rather than a global property of all Galactic 
components. For example, the large negative excursion in the 
40° <£ < 50° bin is due to the W5 1 star formation complex at 
(49.5°, -0.4°). At ^ > 60°, the space density of BGPS sources 
drops sharply, so the effects are dominated by the particular 
distribution of the dense gas within the molecular complex. 
Specifically, in the £ ^ 70° region, there are more dense gas 
structures atb >Q° in the Cygnus X and Cygnus OB 7 region 
than there are at < 0°, while the lower density molecular gas 
does not exhibit this small-number effect. 

7.2. The Flux Density Distribution of BGPS Sources 

In Figure [17] we present two representations of the flux 
density distribution for objects in the BGPS catalog. The 
flux density distribution in 40" apertures is shown in the left 
hand panel and follows a power-law form over nearly three 
orders of magnitude in flux density from the point source 
completeness limit to the upper flux density limit in the sur- 



vey. We have fit a power law to the approximation to the 
differential flux density spectrum dN/dS w AN/ AS oc 
and find that an index a = 2.4 ±0.1 represents the data quite 
well. Despite the similarity to the Salpeter Salpeter ( 1955) in- 
dex a = 2.35, we suspect that this is coincidence because the 
BGPS sources are found at a wide variety of distances. This 
range of distances necessitates a complete treatment of the 
conversion of flux density to mass. Indeed, the Balloon-borne 
Large Aperture Submillimeter Telescope (BLAST) survey of 
the Vela molecular complex at A = 250, 350 and SOQfim found 
a marginally steeper distribution of object mass - dN/dM cx 
M-2^±°-2 for cold cores - than the BGPS finds for flux density 
(iNetterfield et al.ll2009h . A careful treatment of the distances 
to BGPS sources may find the mass distribution to be compa- 
rable between BGPS objects and those found in the BLAST 
survey. 

In the right-panel of Figure [17] we plot the differential flux 
density distribution for the total flux density in the objects S 
(Equation|7]i and find that a power-law is not as good a repre- 
sentation for the catalog. Even though the completeness lim- 
its should be the same for both flux density measurements, the 
power-law behavior for the total flux density only become ap- 
parent above a much higher flux density level of 0.7 Jy where 
a= 1 .9 ± 0. 1 . The difference between these two limits arises 
entirely because of the extended structure accounted for in 
the total flux density measurement. A detection near the com- 
pleteness limit is necessarily similar to a point source. To 
detect extended objects with any degree of fidelity requires a 
much higher peak value of the flux density so that connected, 
extended structure is apparent. Thus, source with flux densi- 
ties near the completeness limit will not sample the full range 
of source sizes that brighter objects will. Hence, the detec- 
tion of extended structure is only apparent in bright sources 
leading to a higher limit. We also note that a power law index 
fit to the distribution above the limit of 0.7 Jy gives a signifi- 
cantly different value compared to the index for flux densities 
derived for 40" apertures. This simply reflects the additional 
effects of the included sizes of the sources. In general, sources 
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with bright peaks (measured in the 40" apertures) also tend to 
be large so their total recovered flux density is larger than for 
smaller sources. This effectively makes the distribution more 
top-heavy, reducing the magnitude of the exponent. 

7.3. Variations in the Flux Density Distribution 

In Figure [18] we show the value derived for the slope of the 
flux density distribution as a function of Galactic latitude over 
the contiguous region of the BGPS. We divide the catalog into 
groups of 50 sources based on the galactic longitude above 
the survey completeness limit of 0.4 Jy. The variation of a 
scatters around the global average of a = 2.4 with marginally 
significant variation as a function of longitude. The variations 
are strongest in regions of high source density suggesting that 
blending may contribute to the derived distribution. Specif- 
ically, the algorithm has a native angular scale over which 
it searches for local maxima (a box with edge length of 2 
beam FWHM). Structures tend to be subdivided on scales 
near this search scale. Large, bright objects in crowded re- 
gions of the plane may be subdivided more often than sources 
in less crowded regions. Indeed, in the I = 60° region, where 
the source density distribution is lowest, the flux density dis- 
tribution is the most shallow. 

7.4. Size Distribution of Sources 

The measured sizes of sources varies significantly for the 
BGPS sources. We plot the distribution of source sizes and 
aspect ratios in Figure [19] for the catalog sources. We show 
the major axis size (?7(T,„a) where 77 = 2.4) determined from 
the moment method and the deconvolved angular radius of 
the source {Or), which includes both major and minor axes 
(Equation|6]l. Because of the deconvolution, the radius 9r can 
be smaller than the beam size. However, we also note that 
some of the major axis sizes are also smaller than the beam 
size. This results from applying moment methods without ac- 
counting for all emission down to the 7 = level. The clip 
at 1(7 truncates the low-significance band of emission around 
the sources. For faint sources, the truncation will clip a sig- 
nificant fraction of the emission resulting in an underestimate 
of the source size. While this effect can be corrected for, we 
refrain from doing so since such corrections require extrapo- 
lating the behavior of the sources below the noise floor in the 
images. 

We also present data for the aspect ratios of the catalog 
sources defined as (Jmaj / crmin- Both the radii and the aspect 
ratios show a peaked distribution with an exponential decline. 
The typical source in the BGPS has a major axis size of 60" 
and a minor axis size of 45". However, there are significant 
numbers of large-aspect-ratio sources which represents asym- 
metric clump structure and, for faint sources, the prevalence 
of filamentary structure in the survey. 

8. SUMMARY 

We present a catalog of 8358 source s extracted from th e 
Bolocam Galactic Plane Survey (BGPS, Ag uirre et ani2009l) . 
The BGPS is a survey of the Galactic Plane from -10° < 
I < 90.5° with extensions into the outer Galaxy. The catalog 
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is generated with an automated source extraction algorithm, 
Bolocat, that is designed for wide-field Bolocam mapping. 
The algorithm mimics a by-eye identification of sources in 
the BGPS. The catalog products are formulated to allow easy 
follow-up since many of the objects in the BGPS are relatively 
unexplored. 

The Bolocat algorithm identifies sources based on their sig- 
nificance with respect to a local estimate of the noise in the 
BGPS maps. Regions with high significance are then subdi- 
vided into individual sources based on the presence of local 
maxima within the region. Each pixel in the BGPS image is 
assigned to at most one catalog source using a seede d water- 
shed, similar to the Clumpfind dWiUiams etal.l l 1994 1) or SEx- 
tractor algorithms (iBertin & Arnoutslll996h . The properties 
of the BGPS sources are measured using the moments of the 
images for each of the assigned sources. 

We have conducted tests of the Bolocat algorithm using ar- 
tificial sources injected into observations in the survey fields. 
The Bolocat algorithm extracts sources with a 99% complete- 
ness limit of 5ct where cr is the local RMS of the noise. Since 
the noise varies across the BGPS fields, the 98% complete- 
ness limit for the survey is 0.4 Jy, set by the highest noise 
fields. 

Using spectroscopic observations of the NH3(1,1) inversion 
transition, we have characterized several of the BGPS sources 
near £ = 32°. We find that the objects in the ca talog are likely 
best described as clumps ( (Williams et al.ll2000 ). though there 
is substantial variation in the properties of the sources. Since 
the catalog is generated using a significance threshold, any 
source that produces a signal > 5(t in the 33" beam size of 
the instrument will be detected. Hence, both nearby, small, 
low-mass objects and distant, large, high-mass objects can be 
detected with similar flux densities. 

We find the flux density distribution of sources in the BGPS 
follows a power-law form. The flux densities of objects ex- 
tracted in a 40" aperture follows a power-law over nearly 
three orders of magnitude: dN/dS^o oc 540 '*^" '. We also find 
that the mean Galactic latitude of objects in the survey lies at 
(b) = -0.095 ± 0.001°, with some variation across the survey 
region according to the presence of large star-forming com- 
plexes. 
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Fig. 17. — Flux density number distribution for objects identified in the BGPS. The left-hand panel shows the flux density distribution for 40" apertures (54o) 
while the right-hand panel shows the total flux density associated with objects (S). We estimate the power law that best fits each distribution for flux densities 
above the peak in the distribution. Note that the source distributions shown are estimators of dN/dS, not dN/d(log M), although the bin widths used to derive 
the distributions vmy logarithmically. 
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Fig. 18. — Power-law index on the flux density distribution as a function of 
Galactic longitude, derived for groups of 50 sources above the survey com- 
pleteness limit of 0.4 Jy. The flux densities are measured in 40" apertures 
across the Galactic plane. 
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Fig. 19. — Distributions of source angular sizes (left) and aspect ratios (right) in the BGPS catalog. In the left panel, the gray shaded histogram indicates the 
major axis size of the sources derived from the moment methods. The black Une-fiUed histogram shows the distribution of the deconvolved radii for the sources. 
The vertical line is plotted at the IJWHM of the beam. In the right panel, the distribution plots the major-to-minor aspect ratio for the moments of the source. 



